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Abstract 

Evaluation of effective or macroscopic coefficients of thermal conductivity under 
coupled heat and moisture transfer is presented. The paper first gives a detailed 
summary on the solution of a simple steady state heat conduction problem with 
an emphasis on various types of boundary conditions applied to the representative 
volume element - a periodic unit cell. Since the results essentially suggest no su- 
periority of any type of boundary conditions, the paper proceeds with the coupled 
nonlinear heat and moisture problem subjecting the selected representative volume 
clement to the prescribed macroscopically uniform heat flux. This allows for a direct 
use of the academic or commercially available codes. Here, the presented results are 
derived with the help of the SIFEL ^ (Slimple Finite Elements) system. 

Key words: Masonry, homogenization, periodic unit cell, coupled heat and 
moisture transfer, effective thermal conductivity 



1 Introduction 

An extensive experimental and numerical analysis of Charles Bridge in Prague has 
been performed only recently to identify the most severe external actions on the 
bridge, see e.g. [17,21,24]. Among others the loading due to spatially and temporarily 
varying temperature and moisture changes appeared to be of paramount importance 
as these effects proved as crucial factors responsible for the nucleation and further 
development of cracks in the bridge. Owing to the complexity of historical masonry 
structures the problem is often addressed in a strictly uncoupled format. Distribution 

^ More information available at http://mech.fsv.cvut. cz/web/?page=software 
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of the temperature field in the macroscopic structure needed for the actual non- 
linear thermo-mechanical analysis can be obtained from an independent solution 
of the macroscopic steady-state heat conduction problem [17]. This step, however, 
requires the knowledge of the macroscopic coefficients of thermal conductivity, and 
possibly their infiuence on actual temperature and moisture content - the principal 
objective of the research presented. 

When dealing with these problems, the application of homogenization techniques 
is inevitable [18]. Solving a set of problem equations on a meso-scale (a composi- 
tion of stone blocks and mortar) provides us with up-scaled macroscopic equations. 
They include a number of effective (macroscopic) transport parameters, which are 
necessary for a detailed analysis of the state of a structure as a whole. A reliable 
methodology of the prediction of these quantities is one of the main goals of our con- 
tribution. Any (multi-scale) approach to coupled heat and moisture transfer draws 
on a cogent description of transport phenomena. An extensive review of this topic 
can be found in [23]. Averaging theories (a micromechanics-based approach), pri- 
marily formulated in [8,9], can be regarded as a counterpart to phenomenological 
ones (a macromechanics-based approach), see [5,2]. Both approaches are explained in 
detail in [16]. Contrary to these trends, phenomenological models are still preferred 
to averaging ones, namely calculating the heat and moisture transfer in building 
materials. 

As outlined in [23], the models for the description of water and water vapor transfer 
can broadly be classified into three main categories, namely convection models, 
diffusion models and hybrid models. The recognized convection model is that of 
Philip and de Vries, see e.g. [19]. A variety of diffusion models were developed 
based on Krischer's original version [12]. A certain drawback of this category lies 
in the absence of cross-effects between the heat and moisture transport. Krischer's 
model was improved by many authors. Let us at least introduce Kiinzel and Kiessl's 
model [14]. Because of the lack of space it is impossible to comment on further 
approaches such as hybrid models or complex models that require the application 
of irreversible thermo- mechanics, [23] . 

While models for transport processes have been developed during several decades, 
the computational methods for multi-scale modeling of these processes in masonry 
on meso and macro scales have emerged only recently. Moreover, most of them are 
prevailingly confined to the effective macroscopic description for heat conduction 
and employ the perturbation method to describe the fiuctuations of temperature 
throughout the heterogeneous material. Different situations are analyzed using the 
homogenization method, which lead to different macroscopic descriptions in [1]. 
Boutin's model analyzing the microstructural influence on heat conduction belongs 
to the same category [4]. It is shown that the higher order terms introduce suc- 
cessive gradients of temperature and tensor characteristics of the microstructure, 
which result in non-local effects. An original approach to homogenization of tran- 
sient heat transfer for some composite materials is proposed in [10]. In that paper. 
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the stochastic second moment perturbation method is used in conjunction with the 
finite element method. Probably the most complex multi-scale analysis for pure heat 
transfer in heterogeneous solids is offered in [18]. The authors established a macro to 
micro transition in terms of the applied boundary conditions and likewise a micro 
to macro transition formulated in the form of consistent averaging relations. See 
also [22] for a similar study with applications to textile composites. 

Homogenization strategies for coupled heat and mass transfer are rather an excep- 
tion, see e.g. [10], and more or less belong to the modeling of a micro to meso rather 
than a meso to macro transition and vice versa. Therein, the condition for a non- 
homogenizable situation, i.e. when it is impossible to find a macroscopic equivalent 
description, is also addressed. 

While setting up the numerical solution of the homogenization problem is relatively 
straightforward, its execution may prove rather complicated if only standard (aca- 
demic or commercial) finite element codes are available. Therefore, apart from the 
determination of the infiuence of temperature and moisture content on effective ther- 
mal conductivities, the present contribution will also be concerned with a strategy of 
predicting the desired effective properties when employing such numerical tools that 
are unable to directly meet the specifics of the homogenization problem formulated 
in the framework of the first-order homogenization approach. To combine the two 
goals of this paper, the following topics will be discussed: 

• In Section 2 we open the subject by reviewing the basic formulas related to first- 
order homogenization clearly identifying the essential similarities between the me- 
chanical and heat conduction problem. Therein and also in Section 3 the steady 
state conditions are assumed in the derivation of individual equations. Certain re- 
strictions on the local temperature field to comply with the assumed macroscopic 
heat fiux or temperature gradient being uniform over a certain representative 
volume element (RVE) are summarized. 

• Numerical treatment of the theoretical formulation is then presented in Section 3 
to provide estimates of the effective thermal conductivities by solving a simple 
steady state heat conduction problem. The main objective of this section is to 
address the effect of boundary conditions on the resulting homogenized properties. 
The fact that these properties are essentially invariant with respect to the choice of 
particular boundary conditions opens the way to the solution of the coupled heat 
and moisture problem using the SIFEL finite element code, a typical representative 
of the finite element software not specifically developed for the solution of the 
homogenization problem. Implementation of the loading and boundary conditions 
into such codes is therefore also discussed. 

• The essentials of theoretical formulation as well as some numerical results are 
presented in Section 4 clearly illustrating the effect of moisture on thermal con- 
ductivity and their dependence on varying material parameters, such as relative 
humidity and initial temperature. Here, unlike the previous sections, the analysis 
is carried out under transient conditions. The presented results suggest that the 
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predicted effective thermal conductivities are almost insensitive to the changes 
in the macroscopic temperature gradient. This is rather appealing, particularly 
in view of the expected use of the homogenized properties in the solution of an 
independent macroscopic steady state heat conduction problem to provide space 
variation of the temperature field within the macroscopic structure (e.g. Charles 
Bridge in Prague), which in turn is used in the actual mechanical analysis [17]. 
• Finally, the essential findings are summarized in Section 5. 

In Section 2, the standard indicial notation is used such that t^i represents a deriva- 
tive of a scalar quantity t with respect to a coordinate Xj. The summation over the 
repeated indexes is assumed. As more typical of the discretized form of the governing 
equations using the finite element method, we adopt in Section 3 the matrix-vector 
notation such that the symbol {a} is reserved for a vector and [A] is employed for 



a matrix representation, respectively [3]. The symbol {Vt} = < — , — , — > then 



represents the gradient of t and {■}"'', [■] stands for the transpose vector or matrix, 
respectively. Finally, the symbol t used in Section 4 represents the time derivative 
of a given quantity t. 



Fig. 1. A representative volume element with macroscopic loading conditions 



2 Fundamentals of homogenization 

Consider a representative volume element in Fig. 1 loaded over the boundary 
r either by the prescribed temperature t derived from a uniform macroscopic tem- 
perature gradient Tj or by a normal heat fiux g^, consistent with the macroscopic 
uniform heat fiux Qi such that 





t = t{x^) +Ti{xi- x°i) onF, 
% = Qi^i on F, 

where z/j represents the i-th component of the outward unit normal. 



(1) 



^ Note that in the case of nonlinear analysis the total quantities would have to be substi- 
tuted by their respective increments. 
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In such a case, the local temperature field admits the following decomposition 



t{x) = t{x'>) + Ti{x^ - xf) + t*{x), (2) 

where Tj are components of the macroscopically uniform temperature gradient vec- 
tor and t*{x) is the fluctuation of the local temperature fleld. Finally, the tem- 
perature t(cc°) at the reference point x^ is introduced to specify the distribution 
of the local temperature field uniquely. This term plays an important role when a 
fully coupled multi-scale analysis is considered and/or when the instantaneous ef- 
fective thermal conductivities depend on the current temperature. If, however, we 
seek only for temperature independent homogenized properties, like in Section 3, 
this term becomes irrelevant and can be safely omitted. 

Eq. (2) immediately follows from the relation between the temperature gradients of 
individual fields 

t,i{x) = T, + t*^{x). (3) 

The micro-temperature gradient averaged over the volume of the representative 
volume element 

yields a scale transition relation, see e.g. [18], 

(t*) = ^ Jj:,{x)dQ{x) = ^ J^r{xMx)dT{x) = 0. (5) 

The boundary integral disappears providing either the fiuctuation part of the tem- 
perature field equals zero (in the case of fully prescribed temperature and/or normal 
heat flux boundary conditions) or the periodic boundary conditions, i.e. the same 
values of t* on opposite sides of a rectangular RVE, are enforced on T. 

When discussing the micro to macro scale transition it is worth mentioning an 
analogy between the basic quantities related to the heat conduction problems (i.e. 
the negative values of temperature gradients (— t,i) or (— Tj) and fluxes qi or Qi as 
their conjugate measures) and the corresponding quantities applied to mechanical 
problems (strains Eij or Eij and the conjugate stress measures (7^ or Ejj). While 
homogenization of mechanical problems benefits from the Hill lemma 

i^ijO'ij) = EijT.ij, (6) 

where Eij, Sj^ are macroscopically uniform strain and stress tensors and e^j, are 
their local counterparts, the Fourier inequality [15] is used to obtain an equivalent 
representation of Eq. (6) given by 

{{-t,i)qi) = -TiQi, (7) 
which is applicable in the case of homogenization of the heat conduction problems. 
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Following the transformation of the volume integral in Eq. (5) to its boundary 
equivalent it is useful to convert the left-hand side of Eq. (7) in a similar way 



t{x)q,{x)dT{x). (8) 



1 

Wt\ Jr 



The term 



-^Jj{x)q,,,{x)dn{x), (9) 

has been omitted in Eq. (8) because of the balance of heat (i.e. qi^i = under steady 
state conditions and the absence of internal heat sources). 

The same boundary conditions, which satisfy Eq. (5), lead to the equivalence of the 
volume averaged microscopic heat flux and the macroscopic heat flux Qt. To prove 
this, the well-known identity p] 

/ qiix)dn{x)= I x,q,{x)dV{x), (10) 

is utilized. The following three situations can be distinguished: 

• In the case of fully prescribed temperature boundary conditions, the subsequent 
combination of Eqs. (2), (8) and (10), in conjunction with the balance condition 

J^q,{x)dT{x) = 0, (11) 

yields 

{t,^q^) = ^ I q.ix) dT{x) + ^ l^ix, - x'^)q,ix) dVix) = 

[ q,{x)dn{x)=T,Qi. (12) 
Jn 

This identity implies 

Q^ = jk\f qi{x)dn{x). (13) 



\n\ Jn 



\n\ Jn 

• A similar situation is met if normal heat flux boundary conditions are prescribed. 
Substituting the prescribed macroscopic flux Qi into Eq. (12) for immediately 
proves the validity of Eq. (13). 

^ f^Xiqu{x)dT{x) = J^Xiqj{x)uj{x)dT{x) = J^{xiqj{x))j dn{x) = J^{5ijqj{x) + 
Xiqj{x)j)d^{x) = f^qi{x)d^{x) 
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• The periodic temperature boundary conditions arc typical of rectangular PUCs. 
From the transition relation (5), the fluctuation temperature t* must be the same 
on the opposite boundaries of PUC. As the anti-periodic normal heat flux applies 
to the periodic boundary conditions, Eq. (12) yields, after certain modifications, 
the expected result (13). 



3 Macroscopic conductivity and resistivity matrices 

Before proceeding with the analysis of a complex coupled nonlinear heat and mois- 
ture conduction problem we review the basic steps for the evaluation of macroscopic 
conductivity and resistivity matrices through the solution of a simple steady state 
heat conduction problem using the finite element method. This allows us to examine 
the influence of various boundary conditions, discussed in the previous section, on 
the predicted homogenized properties and further exploit these results in the next 
section when solving the coupled problem. 

As intimated in the previous section, compare Eqs. (6) - (7), the macroscopic con- 
ductivity and resistivity matrices may be derived by analogy to the homogenized 
stiffness and compliance matrices which apply to mechanical analyses. Therefore, 
we flrst discuss the numerical treatment of Eq. (6), which allows the RVE being 
directly loaded by prescribed macroscopic uniform strain or stress fields Eij or E^, 
resj^ectively. The analogous loading conditions (macroscopic temperature gradient 
Tj or macroscopic heat flux Qi) that apply to heat conduction problem are examined 
next. Such loading conditions are, however, difficult to introduce directly into most 
flnite element codes including the SIFEL program. Instead, the loading conditions 
of the type (1) are needed. This, on the other hand, calls for a special treatment of 
the boundary conditions associated with a fluctuation part of the local temperature 
fleld t*. To elucidate this subject, a brief comment is given in the last part of this 
section. 



3.1 Searching the solution in terms of the fluctuation part of the local temperature 

To begin, consider a rectangular RVE (particular deflnition of an RVE will be given 
later in Section 3.2) subjected either to macroscopic uniform strain {E} or stress {Z}. 
Henceforth, as already mentioned in the introductory part, the standard vector- 
matrix notation is adopted. Limiting our attention to linear elastic theory, thus 
allowing for small displacements, rotations and strains only, provides the virtual 
work representation of Eq. (6) in the form 

{{6ey{a}) = {6Ey{Z}. (14) 
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where the discrctized form of the local strain, in analogy with Eq. (3), and stress 
fields are provided by 

H = {E} + [B.] {ra, W = [L]{e}. (15) 

In Eq. (15) {r*} is the vector of the nodal values of the fluctuation displacement 
field being periodic (the same values of {r*} on opposite sides of an RVE), [By] 
stores derivatives of the displacement shape functions and [L] is the microscopic ma- 
terial stiffness matrix. This equation readily suggests that the solution of Eq. (14) is 
searched in terms of fluctuation rather than actual displacements. Next, introducing 
Eq. (15) into Eq. (14) gives 

{{{6Ey + {6r:y [b,] ^) m m + [b,] {r^})) = {SEyiz}, (le) 

to be satisfied for all kinematically admissible strains {SE} and displacements {5r*}. 

Let {E} be now the prescribed macroscopic strain vector. The above equation then 
reduces to (note that {6E} = {0}) 

([B,]ML][B,]){r:} = -({B,}T[L]){E}, (17) 

to be solved for unknown nodal values of the fluctuation part of the displacement 
field {r*}. The macroscopic constitutive law can be then expressed as 

{Z} = (M) = ([L] ({E} + {B4{r:})) = [L]'^''"^ {E}, (18) 

where [L]'^°"* is the homogenized (macroscopic) material stiffness matrix. Assuming 
e.g. plane-stress conditions the coefficients of the 3x3 matrix [L]'^"'" are defined 
as volume averages of the local fields derived from the solution of three successive 
elasticity problems. To that end, the representative volume element is loaded, in 
turn, by each of the three components of {E}, while the other two vanish. The volume 
stress averages, Eq. (18), normalized with respect to {E} then furnish individual 
columns of [L]'^""*. Note that introduction of the periodic boundary conditions is, in 
this particular case, relatively simple as it is only required to assign, in the FEM 
terminology, the same code numbers to corresponding degrees of freedom of {r*} on 
opposite sides of the RVE. To avoid the rigid body motion and also as a result of 
the application of periodic boundary conditions the rectangular RVE is fixed at all 
corners. 

If, on the other hand, the macroscopic stress {Z} is prescribed (note that in this 
case {6E} ^ {0}) a similar procedure is employed to arrive at a set of the following 
equations for unknown {E} and {r}* 

([L]){E} + ([L] [B,]){r:} = {Z}, 
([L][B,])T{E} + ([B.]T[L][B,]){r:} = {0}. (19) 
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Here, solving for {E} from the three independent elasticity problems and again 
normalizing by {Z} readily provides the macroscopic material compliance matrix 

[M]'^°™ = ([L]'"'-)-i as 

{E} = [M]'"^{T}. (20) 

Derivation of the effective thermal conductivities and resistivities will be now pre- 
sented on the same footing taking advantage of the analogy between mechanical 
and heat conduction problem. To do so, we first write Eq. (7) in the form similar to 
Eq. (14) 

mty{q})^{swTy{Q}. (21) 

The discretized form of the local temperature gradient and heat flux now reads 

{vn = {vr} + [B,]K}, {g} = -[x]{VO, (22) 

where {r*} is now the vector of nodal values of the fluctuation temperature field, [%] 
is the microscopic conductivity matrix and entries of the matrix [Bt] represent the 
derivatives of the temperature shape functions. 

For the prescribed macroscopic temperature gradient {V T} the equality (21) results 
in (compare with Eq. (17)) 

([Bt]"^[x][B,]){r:} = -([B,]-^[x]){Vr}, (23) 

and, successively, in the macroscopic relation for the heat flux 

{Q} = ({?}) = -([X] ({VT} + [B,] {r:})) = - [xr-- {VT}, (24) 

where the homogenized conductivity matrix now follows, with the expected 

restrictions to two-dimensional problem, from the solution of two successive heat 
conduction problems in exactly the same way as already outlined for the homoge- 
nized stiffness matrix [L] The periodic boundary conditions now require the same 
fluctuation temperatures to be ensured on opposite sides of the RVE, see also Fig. 4 
and Eq. (29). In analogy to the mechanical problem a zero value of t* is prescribed 
at all corners of the RVE. 

The prescribed macroscopic heat flux {Q} then provides, in analogy with Eq. (19), 



-(M){VT}-([x] [BM<} = {Q}, 

{[X] [Bt])''{Vr} + ([B,] T [x] [B,]){r:} = {0}, (25) 

{vr} = -[vi/]'"''"{g}, (26) 



where [vl/]'*^"™ = ([x]'^""*) ^ is the effective (macroscopic) resistivity matrix. 
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3.2 Influence of boundary conditions on effective conductivities 

To capture the influence of the selected type of boundary conditions consider an 
RVE displayed in Fig. 2. Such a meso-structure does not represent a typical ma- 
sonry bonding, as in this scheme both the stone blocks and bed and head joints 
are regularly and uniformly distributed. It is also worth noting that such an ar- 
rangement is not perfectly periodic due to the disturbance of periodicity along the 
boundary. This large RVE has been deliberately chosen to demonstrate certain edge 
effects, which may be taken into account when homogenization is carried out using 
commercial and/or academic computer codes. If the RVE was perfectly periodic (i.e. 
one half of the boundary mortar joint was added as sketched by the dashed line in 
Fig. 2), then, due to the infinite periodicity, any stone block along with the adjacent 
part of joints shaded in Fig. 2 would be sufficient for calculating the effective value 
of thermal conductivity. 



? 



I 460 mm 



Fig. 2. Meso-structure consisting of stone blocks and mortar joints 



mortar 


sandstone 


0.87 


1.9 



Table 1 

Phase conductivities xiWm-^K"!] at t = 25°C, (/? = 0% 

In this simple case the coupling effects between the heat and moisture were omitted. 
The phase thermal conductivities are listed in Table 1. First, the effect of imperfect 
periodicity was studied by loading the large RVE in Fig. 2 consisting of 16 stone 
blocks and surrounding bed and head joints by a uniform macroscopic temperature 
gradient V xT = 1 while setting V^T = 0. Fully prescribed boundary tempera- 
tures were considered so that t* = 0. The fluctuation temperature fleld displayed in 
Fig. 3(a) was found from the solution of Eq. (23). Note that the imperfect periodicity 
of the RVE (the edge effect) manifests itself in slightly affecting the periodicity of 
the fluctuation fleld t*. Fig. 3(b) shows similar results derived by assuming the RVE 
in terms of the PUC. Finally, the results plotted in Fig. 3(c) were generated again 
through the solution of the unit cell problem but incorporating the periodic bound- 
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ary conditions. To avoid possible confusion the large RVE, now modified to comply 
with perfect periodicity, was finally examined. Both the infiuence of fully prescribed 
boundary temperatures {t* = on F) and the periodic boundary conditions were 
again studied. 



Large RVE 


Large RVE 


Large RVE 


PUC 


PUC 


non-periodic 


periodic 


periodic 






= on r 


t* = on r 


t*-pcriodic 


= on r 


t*-pcriodic 


1.563 


1.48 


1.47 


1.48 


1.47 


Voight 


1.58 


^hom „m^m i „s^s 

Xx ^ Xx ' ^ Xx 


Reuss 


1.39 


i/x^x'^ = c"'/xf + cyxi 



Table 2 

Effective conductivities Xx°™[Wm~^K~^] 

The corresponding effective thermal conductivities were obtained directly from the 
solution of Eq. (23) together with Eq. (24). Individual values are stored in Table 2 
clearly suggesting, at least in this particular example, the invariance of the predicted 
effective thermal conductivities on the choice of specific boundary conditions. This 
appealing result is further exploited in Section 4.2. The Voight bound derived from 
a simple rule of mixture and the Reuss bound provided by an inverse rule of mixture 
arc presented for illustration. The volume fractions of mortar — 0.31 and stones 
c* — 0.69 correspond to periodic structure (PUC in Fig. 2). 



3.3 Searching the solution in terms of the local temperature - periodic boundary 
conditions in commercial codes 

Recall that the effective thermal conductivities have been obtained assuming the 
solution to be derived in terms of the fluctuation part t* of the local temperature t 
when loading an RVE directly by the prescribed macroscopic temperature gradient 
{Vr} or by the macroscopic uniform heat flux {Q}. The designer, however, must 
often rely on the use of standard flnite element codes, either academic or commercial, 
where the loading is represented in terms of the prescribed boundary temperatures 
or fluxes of the type (1) and the solution is searched directly in terms of the local 
temperatures t{x),x G fl. 

Consider now the most simple case illustrated in Fig. 2, when the RVE is loaded 
by the uniform heat flux \q,,\ = = Qx along the two vertical boundaries at 
X = 0,x = L while no flow boundary conditions are specifled \qiy\ = qy = Qy = 
along the two horizontal boundaries at y = 0,y = H. Eq. (13) is then immediately 
satisfied and the horizontal component of the effective heat conductivity is, in view 
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(b) 



(c) 



Fig. 3. Fluctuation temperature field: (a) large RVE and t* = (b) PUC and t* = (c) 
PUC and t*-periodic 



of Eqs. (26) and (3), provided by 



horn 



(27) 



If no action is taken this result corresponds to the assumption of t* = on F. If, on 
the other hand, the periodic boundary conditions are on demand it is necessary pro- 
ceed as follows. Consider again a two-dimensional rectangular RVE with dimensions 
H and L (see Fig. 4). 



b (stave) 

-Q 



A (master) 



® 



a (slave) 



® 



H 



B (master) 



Fig. 4. Conditions of periodicity 
Observe that for a pair of points (e.g. A - master and a - slave) located on the 
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opposite sides of the PUC the following relations hold: 



tA^ 



( 



dT 
dy 



] 



( 



dT 
dx 



] 



L + 



( 



dT 
dy 



] 



y, + t: + t{x'). 



(28) 



Taking into account the fact that the fluctuation field t* satisfies the periodicity 
condition 



and subtracting corresponding terms on the opposite edges, we finally obtain (com- 
pare with [18]) 



These conditions can be introduced into most commercial software products using 
the multi-point constraint equations. 

Note that neither the assumption of = on F nor the periodic conditions (29) (or 
constraints (30)) are sufficient to specify the local temperature field uniquely. This 

condition is established by introducing a reference temperature t{x^) at an arbitrary 
point 33° G r as indicated in Eqs. (28), recall also the opening Section 2 together 
with Eqs. (l)i and (2). 

In a linear case, when the effective thermal conductivities arc independent of the 
actual temperature, this value is arbitrary and can be set equal to zero in the se- 
lected node of the finite element mesh. In the temperature dependent problem the 
instantaneous effective properties might, however, be strongly influenced by the cur- 
rent temperature. This term then plays an important role. In an example presented 
in Section 4.2 a reference temperature equal to the selected initial temperature was 
assigned to one of the corner nodes of the finite element mesh of the large RVE in 



4 Effect of moisture on thermal conductivity 

So far we restrained our attention to heat conduction problems under steady state 
conditions. To that end, the balance equation, j = 0, Fourier's law and the equiv- 
alence of heat powers on micro- and macro-scales, Eq. (7), were sufficient to derive 
both effective conductivity and resistivity matrices. When analyzing large structural 
systems, such as historical stone bridges, cathedrals and similar historical structures, 
the effect of moisture must be taken into consideration and, therefore, the solution 
of the coupled heat and moisture transfer is desirable. 




(29) 




(30) 



Fig. 2. 
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4-1 Computer codes used for coupled mass and energy transfer 

There exist a number of codes that allow for the description of non-linear and non- 
stationary material behavior as well as cross - effects when studying the moisture 
and heat transfer in heterogeneous infrastructures. Here we recall at least two of 
them. The SIFEL (Simple Finite Elements) computer code developed at our depart- 
ment [13] is a typical representative of programs based on finite element techniques. 
This program was used to derive the results presented in Section 4.2. The DELPHIN 
program [6,7] is on the other hand developed on the bases of finite control volume 
method (FCVM). Although not specifically used in this contribution this program 
has proved useful in a number of verification tests. 

In this section we mention two particular models implemented in the SIFEL pro- 
gram. The model of Lewis and Schrefier is likely the most advanced hygro-thermo- 
mechanical model so far employed for practical calculations of transport processes 
in soils and other deformed porous media. As its detailed description can be found 
e.g. in [16], we just summarize the basic features of this model. The complete set 
of model equations comprises the linear balance (equilibrium) equation formulated 
for a multiphase body, the energy balance equation and the equations of mass con- 
servation for liquid water and gas. Assuming the rigidity of the solid phase, the 
number of fundamental unknowns in the theory of the coupled moisture and heat 
transfer reduces to three, namely the temperature t and the pressures in the hquid 
and gaseous phases p^,p^, respectively. These unknown functions can be solved from 
three basic equations 

• one energy balance equation, 

• two equations of mass conservation. 

The second model implemented in the SIFEL code is that of Kiinzel and Kiessl [14]. 
It is much simpler then the previous one and, despite the fact that the cross-effects 
are missing in this model, it describes all substantial phenomena and its results 
comply well with experimentally obtained data. For these reasons it has been chosen 
to carry out the case study on masonry microstructure and will be commented on 
in more detail in this paragraph. 

Vapor diffusion and liquid transport due to differences in capillary suction stress 
are the main mechanisms of moisture transport. Although an interaction of liquid 
and vapor fiuxes under very humid conditions cannot be excluded, they are treated 
as independent processes in this model (vapor diffusion is most important in large 
pores, whereas liquid transport takes place on pore surfaces, in crevices and small 
capillaries) . 

Thermo-diffusion is neglected in this model and vapor flux, {^v}, is described by the 
following equations 

{q.} = -Sp{Vp}, (31) 
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where Sp is the vapor permeabihty of the porous material and p is the vapor pressure. 
The hquid transport incorporates the hquid flow in the absorbed layer (surface dif- 
fusion) and in the water filled capillaries (capillary transport). The driving potential 
in both cases is capillary pressure (suction stress) and/or, according to Kelvin's law, 
relative humidity (f. Hence, the flux of liquid water, {^w}, can be written as 

{gw} = -D^i^ip}, (32) 

where I?^ = I?^-^ is the liquid conductivity , D^, — Dyj{w/wf) is the liquid 

diffusivity , dw / dy? is the derivative of the water retention function and w/ Wf is 
the water content related to the capillary saturation. Finally, the heat flux, {^t}, 
attains the following form 

{ft} = -x{Vt} - /i„5p{V(mat)}, (33) 

where hy is the speciflc heat of evaporation and Psat is the saturation vapor pressure. 

In this model, the resulting set of differential equations for the description of heat 
and moisture transfer, expressed in terms of temperature and relative humidity ip, 
assumed the following form: 

• The energy balance equation 

[pc + ^) i = {VY \x{Vt} + /i„5p{V((^Psat)}] (34) 

• The conservation of mass equation 

= {Wy [D^{W^} + 5p{V((^p3at)}] . (35) 

In Eqs. (34) - (35) p is the material density, c is the specific heat capacity and 
Hr^j is the enthalpy of material moisture. Details concerning the specification of 
material parameters can be found in [14]. The storage terms appear on the left 
hand side of Eqs. (34) and (35). The conductivity heat flux and the enthalpy flux 
by vapor diffusion with phase changes in Eq. (34) are strongly influenced by the 
moisture fields. The vapor fiux in Eq. (35) is governed by both the temperature and 
moisture fields due to exponential changes of the saturation vapor pressure, Psati 
with temperature. Vapor and liquid convection caused by total pressure differences 
or gravity forces as well as the changes of enthalpy by liquid flow are neglected in 
this model. 

In the SIFEL program these equations are converted to a weak formulation. Using 
the weighted residual statement the energy balance equation becomes 
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Fig. 5. A two dimensional region (e.g. an RVE in Fig. 2 with prescribed heat flux) 



c+ ^] i - {vy [x{vt} + KS,{v{ipps,,)}]\ dn+ 



+ I wt {{uy [x{Vt} + K5p{V{cpps,,)}] + q,,) dr = 0, (36) 

where Wt is the weighting function such that = on F^, see Fig. 5. Applying 
Green's theorem then yields 



C + i + {Ww,y [X{V0 + /i.5p{V((^Psat)}]| df]+ 

+ [^Wtqt^dT^O, (37) 

where is the prescribed heat flux perpendicular to the boundary Ff. On Fj the 
temperature t is equal to its prescribed value t. Vector {z/} stores the components 
of the unit outward normal defined already in Section 2, see also Fig. 5. 

The weak formulation for the conservation of mass can be derived analogously 



In {^'^ ~ ^^^^ [DA^'p} + <^p{v(mat)}]^ dn+ 

+ {{uy [D^{V^} + 5p{V(^Psat)}] + q^.) dF, (38) 

where w^p is the corresponding weighting function such that = on F<^. Appli- 
cation of Green's theorem finally gives 



+ / , ^f^wu dr, (39) 



where q^^, is the prescribed flux of water perpendicular to the boundary F^. On F^ 
the relative humidity is again equal to its prescribed value (p. 
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A detailed specification of boundary conditions used in our particular case will be 
given in the next section. For details on numerical integration of Eqs. (37) and (39) 
the interested reader is referred to [13]. 



4-2 Effective conductivities as a function of relative humidity and initial tempera- 
ture 



As an example the non-periodic 16-blocks RVE in Fig. 2 was examined to assess 
the influence of moisture and initial temperature on the effective thermal conduc- 
tivities. Transient heat conduction problem was solved with the help of SIFEL finite 
element program exploiting the Kiinzel and Kiessl model from the previous section. 
The loading and initial conditions are plotted in Fig. 6(a). Although there are no 
restrictions to introduce the periodic boundary conditions in the SIFEL program 
through Eqs. (30), the present analysis assumed, with reference to Section 3.2, the 
fiuctuation part of temperature t* = on F. This simplification on the other hand 
allows for a direct comparative study using other codes such as DELPHIN, where 
implementation of periodic boundary conditions might be precluded. 
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Fig. 6. Example: (a) Geometry and boundary conditions, (b) Evolution of temperature at 





X 



Parameter 


Mortar 


Sandstone 


density [kgm~^] 


1700 


1964 


specific heat c [Jkg~-'^K~^] 


1000 


900 


vapor diffusion resistance number /x, see Eq. (44) 


12 


10 



Table 3 

Constant material parameters of individual phases 



The constant material parameters are listed in Table 3. The moisture dependent 
variation of phase thermal conductivities and liquid diffusivities is plotted in Fig. 7 
together with sorption isotherms obtained from the following approximations [7] 
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w 



a/1 - 'Phyg' 

, ^ - Vhyg , V 

W = Whyg + [Wsat - Whyg) , 

J- - (Phyg 



(40) 
(41) 



where the saturation moisture content Wgat = 0.3, the hygroscopic moisture content 
Whyg = 0.02 and the hygroscopic relative humidity iphyg = 95% were assumed. 
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Fig. 7. Example - material data: (a) variation of phase thermal conductivities as a function 
of volume moisture, (b) variation of phase liquid diffusivities as a function of volume 
moisture, (c) phase sorption isotherms 



The remaining parameters that enter Eqs. (31) - (35) are provided by 
• hv specific heat of evaporation [Jkg~^] 



K = 2.5008 • 10^ 



273 15 \ (0-167+t-3.67.10-4) 

~t ) 



(42) 



• Psat saturation vapor pressure [Pa], see [11] 

4042.9 



Psat = exp (^23.5771 



t - 37.58 



(43) 
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5p vapor permeability of the porous material [kgm ^Pa 



(44) 



where yU is the vapor diffusion resistance number and 5 is the vapor diffusion 
coefficient in air [kgm^^s^^Pa^^] given by, see [20] 



5 



2.306 ■ 10-5 ■ J9a / t 



R^-t-p V273.15 



1.81 



(45) 



with p set equal to atmospheric pressure pa = 101325 Pa and = R/M^ = 
461.5; R is the gas constant (8314.41 Jkmol-^K"^) and My^ is the molar mass of 
water (18.01528 kgmol'^). 

The total temperature field t = T+t*, which corresponds to = lOOWm"^, appears 
in Fig. 8(a). In this case the room initial temperature tj„ = 25° C and relative 
humidity cp = 60% were considered, recall Fig. 6(a). As evident from Fig. 8(b), 
the volume moisture field w (m^/m^) considerably varies throughout the RVE thus 
affecting the macroscopic thermal conductivity slightly depending on the size of the 
RVE. 




(a) (b) 

Fig. 8. Example: (a) Distribution of temperature (b) Distribution of volume moisture 

Several interesting results have been derived within the scope of our numerical ex- 
periments. The graph in Fig. 9(a) displays the dependence of the macroscopic heat 
conductivity -^^""^ on the average relative humidity (p. Individual values of ■^'^""^ 
follow again from Eq. (27). Note that this equation corresponds to steady state con- 
ditions attained after running the transient heat conduction problem for about 50 
hours, see Fig. 6(b) showing an evolution of temperature at a reference point 
with respect to time. It should be pointed out that the nearly linear dependence 
between these two variables at a sub-hygroscopic region (up to about 95% - 98%) 
becomes highly non-linear once the hygroscopic moisture is attained and approaches 
the capillary and even maximum (vacuum) saturation. Fig. 9(b) displays the rela- 
tion between the macroscopic thermal conductivity and the initial temperature for 
different levels of the relative humidity (fin- 

These information were subsequently utilized when performing the analysis of Charles 
Bridge. It has been observed experimentally that different sections of the bridge ex- 



19 



perience different moisture content. The results in Fig. 9 thus allowed us to assign 
to each measured value the corresponding value of in the macroscopic steady 
state analysis of Charles Bridge. This then provided the space distribution of tem- 
perature eventually used in the thermo-mechanical analysis of the bridge [17]. 
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Fig. 9. Example - evolution of x'"'™ as a function of: (a) cp at tin = 25°C, (b) Un 

The sensitivity of effective thermal conductivities on the variation of macroscopic 
temperature gradient was also studied. Clearly, the results presented in Fig. 10 sug- 
gest that the coefficient of the macroscopic thermal conductivity -^^""^ is almost 
insensitive to the changes in the macroscopic temperature gradient Tj, which con- 
siderably simplifies the micro (meso)-macro approach. 
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Fig. 10. Example: sensitivity of x'*""* on variation of Tj 



5 Conclusions 



Homogenization techniques applied to a meso-scale provide the values of effective 
transport parameters, such as the coefficient of thermal conductivity (and/or mois- 
ture permeability), and allow for their dependence on both macroscopic temperature 
and moisture fields. As efficient tools, the SIFEL and DELPHIN computer codes can 
be effectively used when analyzing transport processes within a certain PUC to get 
distributions of these fields on a meso-scale, while taking into account the coupling 
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Colour image Binary image SEPUC 




Fig. 11. Construction of statistically equivalent periodic unit cell (SEPUC) for random- 
coursed masonry [25] 



effects in transport phenomena. With reference to thermal conductivity the following 
conclusions can be pointed out: 

• At low levels of moisture the coefficient of thermal conductivity varies almost 
linearly with the relative humidity ^9. This relation becomes strongly non-linear 
as soon as the hygroscopic moisture (9? ~ 95%) is attained. 

• The macroscopic volume moisture (m^/m^) extensively varies throughout the RVE 
thus affecting the value of macroscopic thermal conductivity. 

• There exists a certain dependence of the macroscopic thermal conductivity on the 
initial temperature tin- But it is distinctive just at higher temperatures, say above 
30°C. 

• A very important finding is that the macroscopic conductivity is nearly insensitive 
to the macroscopic temperature gradient, even in the non-linear range, which 
simplifies the meso-macro approach. 

Finally recall a relatively good correspondence between the numerically obtained 
values of the effective conductivities and those calculated from simple rules of mix- 
ture, see Table 2. It should be, however, mentioned that such a good agreement is 
merely attributed to the selected geometrical arrangement of the two phases. When 
more realistic periodic unit cells (such as the one displayed in Fig. 11) were used, 
the simple bounds were too far apart to be of any practical use. 
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